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#-/-\ ' The effect of geometrical shape of eroding absolutely rigid particles on the threshold rate of 

f— ^ , failure has been studied. The Shtaerman-Kilchevsky theory of quasi-static blunt impact, which 

£SJ ' generalizes Hertzs classical impact theory, is used for modeling the frictionless contact interaction 

of an axially-symmetric particle with an elastic half-space. The incubation time fracture criterion 
is applied for predicting surface fracture. It is shown that there exist a critical value of the particle 
shape parameter such that for all its lower values the fracture energy possesses a nonzero minimal 
J^ , value. 

5_i , 
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1 Introduction 

Surface fracture of materials due to solid particle erosion is encountered in a wide range of applications, 
including aerospace applications and degradation of components used in machinery related to mining 
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and oil-drilling [1 ] . In recent years a number of experimental [2J [3J and theoretical [H [51 [5J studies on 
the behavior of hard surfaces when subjected to solid particle erosion with high-impact velocities have 
been published. However, theoretical modeling of erosion phenomenon has still to be developed to get a 
deeper insight into the brittle fracture mechanisms of erosion at high impact velocities. 

In analyzing dynamic strength of materials in a wide range of loading rates from quasi-static to 
extremely dynamic, one often encounters a contradiction between available experimental data and the 
classical approach for fracture characterization that is based on quasi-static consideration of the failure 
process 0151 |S]. In particular, a number of well-established experiments demonstrate that under high- 
rate dynamic loads, different materials are able to sustain loads substantially higher than those that 
cause fracture initiation in quasi-static conditions [5J . Study of this phenomenon led to the development 
of a concept of incubation time [TU1 [5J . 

Most studies dealing with problems of erosion fracture jll[ 112] develop a phcnomcnological approach 
and usually employ fracture criteria based on extrapolation of quasi-static fracture criteria to dynamic 
conditions. However, the temporal characteristics of the fracture process are usually not taken into con- 
sideration, though sometimes even inertia and temporal characteristics of the applied load are accounted 
for. In the present study, an incubation time fracture criterion [TU] is used for predicting surface fracture 
in the erosion process. Another assumption that is usually made is that the eroding particles are con- 
sidered to be spherical. This simplifying assumption is almost always violated in reality (in particular 
for sand particles). The objective of this work was to develop a mathematical modeling framework for 
studying solid particle erosion in a wide range of particle impact velocities with varying particles shapes. 
One of the main results is that there exist a critical value A* for the particle shape parameter A such 
that for all lower values of A the fracture energy possesses its nonzero minimal value. 

The paper is organized as follows. In Section [2] we consider the quasi-static theory of blunt impact 
as it was developed by Kilchevsky [13] based on the theory of axisymmetric frictionless contact for a 
polynomial indenter due to Shtaerman [Mj. We generalize the Shtacrman- Kilchevsky theory for the case 
of a general non-polynomial axisymmetric indenter according to the solution obtained by Galin |15] . The 
final form of force-displacement relationship for the blunt indenter was obtained using recently published 
results by Borodich and Keer [16]. In Section [3l applying the approach developed by Johnson [17], we 
consider the question of indentation fracture beneath a blunt indenter. The geometry of eroding particles 
is described in Section 01 An incubation time-based fracture criterion [TU] is introduced in Section [5J 
Prediction of the threshold fracture energy in the erosion with equiaxed supcrellipsoid particles is given 
in Section [SJ The developed mathematical model is illustrated in Section [7] for the set of material 
properties parameters taken from the literature [BJ . A discussion of the results obtained has been also 
outlined in Section [7] Finally, in Section [8j we formulate our conclusion. 

2 Shtaerman — Kilchevsky theory of quasi-static blunt impact 

Let us consider quasi-static indentation of an clastic half-space by a blunt indenter with the shape 
function 

z = Ar x , (1) 

where 1 < A is a real number, A is a constant having the dimension [L 1_A ] with L being the dimension 
of length, r is a polar radius. 

According to [T5] (see also [IB]), the relationship between the contact force, P, and the indenter 
displacement, w, can be represented as 

P = k l w h i k , (2) 

where 



Here, E and v being Young's modulus and Poisson's ratio of the elastic semi-infinite body, T{x) is the 
Gamma function. 

In the case A = 2n, where n is a natural number, according to |14| . we will have 

E j. An f (2n-l)l\ \^ 

kx -T~^ A 2 "2^TT^ (2n)!! ) ' (4j 

where (2ra - 1)!! = 1 ■ 3 . . . (2n - 1) and (2n)!! = 2 ■ 4 . . . {In). 

Let m denote the mass of the indenter. Then, in view of ([2]). the equation of the indenter motion 
can be written as follows: 

m— j- = -fat/A (5) 

Here we introduced the short notation 

Equation ((SJ) should be supplied with the initial conditions 

«;(0) = 0, ^(0)=t*, (7) 

where vq is an initial velocity of the indenter. 
The energy conservation law yields 

j(v 2 -v 2 o) = -£- 1 ™ 0+l , (8) 

where v = dw/dt is the indenter velocity, and the initial conditions ([7]) were taken into account. 
Putting v — in Eq. (|8|) . we get the maximum penetration 

The impact duration To = 2to, where to is the time instant at which the penetration w{t) reaches its 
maximum wq. From Eq. (|8|), it follows that 

to = —b, (10) 

vo 



7 <3 = / n = , n . nw 0+3 \ ' y 11 ' 



where 

V/3+ 

/o ^^^^TT G9 + 1)r (_gt3 y ) 

Note that formulas ((9j) and ([TT]) were also derived in [18]. 
Let us introduce the following dimcnsionless variables: 

u = — , T = 1. (12) 

w w 

Making use of the new variables (|12[) . we rewrite Eq. ([5]) and the initial conditions (JT)) as follows: 

£ + <2^.c da, 

«(0) = 0, f(0) = L (14) 
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The unique solution of Eq. (|T5|) satisfying the initial conditions (|14l) will be denoted by U(j3; r) 
Thus, the dependence on time of the penetration w(t) is given by 



(t) = u; Cr(/3; — t), (15) 



where te [0,2*o]. 

Integrating once Eq. (|T3]) and taking into account the initial conditions (fl4)l . we get 



^V+^ +1 = l. (16) 

In the penetration stage, when u(t) increases, Eq. (|16[) yields the following equation: 



o VT^W 1 



= t. (17) 



The above relationship represents an implicit equation for u(t). Indeed, let F(u) denote the left-hand 
side of Eq. (fTT)) . Then, u(t) = F _1 (t), where i* 1-1 is the inverse function of F. Note also that Ip = F(l) 
(see Eq. (HI])). 

3 Indentation fracture beneath a blunt indenter 

The contact pressure beneath a blunt indenter with the shape function |T]) is given by 



P(r) = i -^ 1 P, / '""V +£ 2 )^C (18) 



2 jo 

Here, p = r/a is the dimcnsionless radial coordinate, po = P/{~Ka~) is the mean contact pressure, a is 
the radius of the contact area given by 

/ (i-^)P(A + i)r(A) \^ 

In the case A = 2n, according to |14j . the contact pressure can be represented in the integral form as 
follows: 

(2n + l) 



where S n (p) is the n-th Shtaerman polynomial defined as 



(2n 2)!! / 2n-2 , ! „2n-4 , 3 n 2«-6 

Mp) " (2n-l)!!^ + 2 P + 2~4 P 

(2n-5)!! 2 (2n-3)!! 



(2n-4)!r (2n-2)!! / 
Correspondingly, the contact radius a is related to the contact force P by the equation 

(l-v 2 )P(2n+iy.\\^+~ 1 
EA 4n(2n)! 



(21) 



(22) 



The radial and circular stresses throughout the surface of the clastic half-space can be expressed in 
terms of the function 



as follows \T7\ 
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According to Eqs. (|24p and (|2"5]) . at the edge of the contact area, we will have 



K Ja 



(26) 



Now, substituting the expression (|23|) into the right-hand side of Eq. ([26)) and changing the order of 
integration, we obtain 

(27) 



(28) 



a r {l) = (l-2v) / r)P(v)dv- 
Jo 

Finally, by the definition of the contact force, we get 

(1 - 2v)P 



<r r (l) 



27ra 2 



Note also that crg(l) < 0, or, in other words, the circular stress at the edge of the contact area is 
compressive. The stress field in an clastic half-space beneath an arbitrary axisymmctric indenter was 
studied in detail in [15] . 



4 Geometry of an axisymmetrical superellipsoid particle 



A general axisymmetrical superellipsoid centered in the coordinate origin is described by the following 
implicit equation: 

/ / o . ~\ A 

Z 



\j x 2 + ;" 2 x 



//- 



B 



C 



= 1. 



(29) 



The shape of the superellipsoid (|29|) around its pole will be described by Eq. (p} up to the terms of 
higher order, if the following relation takes place between the constants A, B, and C: 



A = 



C 
XB^' 



(30) 



With A and A being fixed constants, the shape of the supcrellipsoid ([2T))) will be defined uniquely, if 
we require additionally that the volume, V, of the supcrellipsoid is also fixed, that is 



47r„ „rmr(± 

'3 

.A 

where T(x) is the Gamma function. 



v = ^-b 2 c va ;,; a; , (31) 

3A r(^ 



5 Incubation time-based fracture criterion 

In this section, an incubation time fracture criterion is used to predict fracture in the case of impact 
interaction of a superellipsoid particle with an elastic half-space. In the simplest case, the incubation 
time fracture criterion is formulated as follows [ID] : 

1 r* 

a r (t') dt' = a c . (32) 

Here, <j c is the tensile strength of the elastic material (this parameter is evaluated for quasi-static loading 
conditions), r is the incubation time of the fracture process. 

Let uj(t), t £ [0,To], be the shape function of an impact loading pulse. Then the stress variation is 
given by 

er,.(i) = a max w(t), (33) 

where er max is the stress amplitude. 

Substituting the expression Q33p into Eq. (|32[) , we obtain the critical (threshold) amplitude, c^ax; 
leading to fracture, i.e., 

. (34) 



1 
Lo(t') dt' 



te[o,T ] 



t—T 



Here, Tq is the duration of the impact process. 

According to Eqs. ©, ©, (EH]), and ([19]). we will have 



a r {t) = H 1 (A) (1 1 2 "l E Aiw(t)^, (35) 

1 — v* 



where we introduced the notation 



Hi (A) = ^rrv ; ^r- (36) 

iV ' 2tt A + l r(A)i 



On the other hand, in view of Eqs. ([21), ©, and ([T5J) . we get 



w(t)=w U[^±±;^t), 

A wo ' 



WQ = Tl 2 (\)(mv 2 )^A^TT(^-y* + \ (38) 



where 

'(2A + i)(A + i)\^ FT /r(|) 2, 



^W = [ ^ ' ) [^ I (39) 



Now, substituting (|57|) and (j3"8")l into Eq. (|35p . we obtain 



<W = n 3 (A)(l - 2^)( T — ^ 2A+1 ^^Vr( mw 2) wr ; (40) 

V A w J 



A-l 



113(A) = IIi (A)II 2 (A)-k-. (42) 

In view of the notation ©, Eq. (|4"Tj) can be rewritten as follows: 

u{t) = U*-?(p-At). (43) 

Since the shape function (|43[) of the impact loading pulse is symmetric with respect to the line t = to, 
where to = To/2, Eq. (|M)) can be transformed into the following one: 

er = (44) 

max ft +T/2 ' *> ' 

/ Uj{t') dt' 

Jto-r/2 

Let us now introduce the notation 

T(/3,i )= [ t0+T/2 U^((3;^t')dt'. (45) 

Jto-r/2 V Wo / 

The evaluation of the integral (|4"5j) depends on the relation between the impact duration To = 2to 
and the incubation time. Indeed, if 2to < t, then 

T(/3, to) = 2 / ° t/ 2 -' 3 (>; ^i') dt' . (46) 

Jo V w / 

By the change of the integration variable to s = vot' /wo, we simplify the integral ([4"fJ|) as follows: 

T(/Mo) = ^ [ l> U 2 -P(p;s)ds. (47) 

J/3 Jo 



Now, if r < 2to, then in the same way we obtain 

T(/3,to) = ^r C/^^s)^. (48) 

Recall that the integral Ip is defined by formula (jTTJ). 



6 Prediction of the threshold fracture energy 

In what follows, we consider the case of an equiaxed supcrcllipsoid particle. In other words, along with 
the assumptions (pRJ)) and (|3"Tj) . we put 

C = B. (49) 



Thus, according to Eqs. ([30]) . (|3Tj) . and (|49|) . the particle's mass is given by 

m = n 4 (A)p ^ _ ^ T , (50) 
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where po is a parameter of load intensity, having the dimension of mass density. 

n lW = f^r(f)r(I)r(l)-. (51 , 

In view of ([50]) and (|40p , the maximum tensile stress generated in the impacted medium is evaluated 
as follows: 

<W = n 5 (A)(l - 2u)( Y ^)^ T pf^vf^. (52) 

Here we introduced the notation 

n 5 (A) = n 3 (A)n4(A)^TT. (53) 

Observe that the right-hand side of Eq. (152"]) docs not depend on the geometrical parameter A. 
Further, according to Eqs. (|TD]), ©, and (J5U)) . we obtain the penetration duration to in the form 

t = n 6 (\)(±^)^p^A-^v-^-\ (54) 



where 

n 6 (A) = ^n 2 (A)n 4 (A)^TT. (55) 

Following [6], we extract the quantities A and vq from the system (|52|) . (|54]l. First of all, Eq. f|52[) 
immediately yields 

" 0= ( v ( i^2,)n 5( A) J l-B-J Po "' (56) 



Second, substituting the expression ([56]) for «o into the right-hand side of Eq. (|54|) . we obtain after 
some algebra the following equation: 



A *■-! = 



n 6 (A)V(l-2z,)n 5 (A)J U-W ^0 • (57) 



Evaluating the initial kinetic energy, eq, of the supcrcllipsoid particle, one can estimate the energy 
required to create fracture in the impacted medium as follows: 

o 4X+5 

_ fg(<w) '(*-*> 

£ — "A 3 10 _ A • (,00) 

Here, a\ is a dimensionlcss coefficient given by 



10-A 



H 4 (A) (i-^)^iy 

^ ITT l" \ V3 4A + 5 ' V 0y / 

211 6 (A) ((i_2^)n 5 (A)) 3 <^ T ) 
Note that a\ depends on Poisson's ratio (this is not indicated in notation for simplicity). 

7 Numerical results and discussions 

As shown in a number publications [20l [TOJ [21] [22j [23] , the criterion (|32|) can be successfully used for 
predicting the fracture initiation and structural transformations in brittle solids under dynamic loading. 
For slow loading rates (and, correspondingly, times to fracture that are essentially bigger than t), 
the condition (|32p for crack initiation gives the same predictions as the classic critical stress criterion. 



For high loading rates (and, correspondingly, times to fracture comparable with r), all the variety of 
effects experimentally observed in dynamic fracture experiments (9[ [20] can be explained using (1321) , 
both qualitatively and quantitatively. It should be emphasized that the criterion itself gives correct 
predictions in a wide range of loading rates, and it is not necessary to pay a special attention to the 
particular time scale of the impact problem. 

In our numerical calculations, the properties of the impacted material are taken to be equal to the 
properties of gabbro-diabase [6] as follows: E = 6.2 x 10 9 N/m 2 , v = 0.26, er c = 44.04 x 10 6 N/m 2 , 
r = 44 [is. The mass density of eroding particles is po = 2400 kg/m 3 . 
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Figure 1: The initial kinetic energy, Eq, as a function of the impulse length, 2i , for different shapes of 
the blunt indenter determined by the value of the parameter A indicated in the legend sections (Fig. [TJt 
and Fig.QJ>). For all A < 5.5 the energy possesses its nonzero minimal value. In the critical case A = 5.5 
the minimal energy corresponds to the horizontal part of the graph. 



First of all, observe (see Fig. QJi) that with decreasing the shape parameter A, the fracture energy £ 
decreases too. This dependency is so strong that a logarithmic scale is required to show the trend. Recall 
that the case A = 2 corresponds to spherical particles, while with decreasing A, the contact surface shape 
approaches a cone (with the limit value A = 1). It is clear that infinitesimally small fracture energies 
are required for a sharp indenter, which produces a singular contact pressure distribution. By the way, 
the same can be said of a cylindrical indenter with a sharp edge. Note also that under the assumption 
of elastic impact the fracture energy coincides with the kinetic energy of eroding particles. 

Remark. It is interesting to observe that for all A < 5.5 the fracture energy possesses a nonzero 
minimal value. In other words, there exists an energetically optimal mode of impact erosion such that 
the energy input for fracture is minimal. From the practical point of view, this effect can be used, for 
example, for optimization of the polishing process so that the energy cost of the process is minimized 
[2"4"] . After determining the critical value £q, for which the fracture energy eo takes its minimum, Eq. (|54p 
allows one to determine the corresponding value Vq for the eroding flow velocity. In view of Eq. (|50[) . we 
will have 



t* 



^n 2 (A) 



i 



E 



m 2TFIyl2Hl( Wo *) 2J+T, 



Because t$ depends on A, the optimal value Vq will also depend on this shape parameter as well as on the 
material parameters E and v, mass of eroding particles m and on the geometrical parameter A. In view 
of Eqs. (|30|) and (|49|) . we have the relation A = A - 1 C 1 , where C is the particle semi-axis. It should be 
noticed, of course, that the eroding particles do usually have a wide range of size and shape parameters 
and this circumstance should be taken into account. Also, friction docs play a role, especially, for oblique 
impacts [25) . In the case of normal impact incidence, as a first approximation, friction effects can be 
usually neglected. 



However, it is more interesting to find that there is the so-called critical value A* = 5.5 of the shape 
parameter, when the minimal energy corresponds to the horizontal part of the graph (see Fig. [T]d). This 
means that for A = A* , the fracture energy remains the same for small values of the duration of loading 
impulse 2tg comparable with the incubation time r. 
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Figure 2: The initial kinetic energy, eq, as a function of the impulse length, 2io, for the shapes of the 
blunt indenter determined by the parameter A larger than its critical value A*, that is A > 5.5. There is 
no minimum for the energy (the shorter impulse the lower energy level is accessible). 

Fig. [2] illustrates the influence of the impact duration 2io on the fracture energy £q for blunt particles 
with the super-critical shape of the contacting surface. The absence of the minimum can be explained by 
the fact that the contacting surface with increasing its bluntness leads to the contact stress concentration. 
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a) 
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Figure 3: Fig. [3ja): The minimum of the initial kinetic energy, e$, in the logarithmic scale as a function 
of the parameter A determining the shape of the blunt indenter. Fig. [3jb): The length of the impulse, 
f"minj corresponding to the minimal value of the initial kinetic energy, e m i n . 

Finally, Figs. [3^, and [3}} present the variation of the minimum £ m m of the fracture energy eq and 
the corresponding impact duration r m i n with variation of the shape parameter A. It is seen that with A 
approaching its critical value A*, the optimal impact duration r m ; n approaches the incubation time r. 
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8 Conclusion 

In the present paper, the effect of geometrical shape of eroding absolutely rigid particles on the threshold 
rate of failure has been studied. The eroding particles are assumed in the form of axisymmctrical cquiaxed 
supcrcllipsoid, and the case of normal impact is considered. The dynamic contact interaction between 
an axially-symmetric particle with an elastic half-space has been described in the framework of the 
Shtaerman - Kilchcvsky theory of quasi-static blunt impact. Stress field created at the contact surface 
as a result of the impact interaction is estimated using Galin's solution of the contact problem for a blunt 
indentcr. For predicting surface fracture beneath a frictionless blunt indenter, we employ the incubation 
time fracture criterion. 

The effect of geometric shape of eroding particles on the threshold rate of failure was very recently 
considered in [35] , where the following three special cases were studied: spherical [2U |B] , cylindrical [57] , 
and paraboloidal shape with A = 4. In all papers mentioned above, the particle motion function was 
assumed according to the sine law, while in the present analysis the particle equation of motion is solved 
implicitly. However, the main advantage of the employed approach concerns the developed mathematical 
modeling framework for studying the geometry effect in the whole range of the particle shape parameter 
A. In the case of spherical particles, the present analysis reduces to the theory developed in [2U [5] with 
the only difference that the actual particle motion slightly differs from the sine law. However, in the 
other two cases, the developed approach differs also by the use of the cquiaxed supcrcllipsoid shape for 
eroding particles. 

Based on the present analysis, the kinetic energy necessary for a rigid particle to fracture the impacted 
surface is analyzed. Threshold particle velocity (and, hence, threshold kinetic energy) corresponding to 
the initiation moment of surface fracture in the elastic half-space is determined. As the main results of 
the present study it is shown that the value of the fracture energy does significantly depend on the load 
duration and has a marked minimum in the so-called subcritical case when A < A* . It is also shown that 
A > A* the fracture energy achives its zero minimal value with decreasing impact duration. Existence 
of energetically optimal modes of dynamic impact is claimed. In particular, it is demonstrated that the 
energy input for fracture can be optimized so that the energy cost of the process is minimized. The 
effect itself turns out to be dependent on the bluntness of the particle's contacting surface. 

Acknowledgements 

LA. and Yu.P. were supported by the Russian Foundation for Basic Research (project Nos. 10-08-00966- 
a and 1 1-01-0049 1-a); G.M. gratefully acknowledges the support from the European Union Seventh 
Framework Programme Marie Curie Programme (project No. PIAP-GA-284544-PARM-2). 

References 

[1] B. Arcfi, A. Scttari and P. Angman, Wear 259 (2005) p.263. 

[2] D.W. Wheeler and R.J.K. Wood, Phil. Mag. 89 (2009) p.285. 

[3] A. Gnanavelu, N. Kapur, A. Neville and J.F. Florcs, Wear 267 (2009) p. 1935. 

[4] J.M. Fan, C.Y. Wang and J. Wang, Wear 266 (2009) p.968. 

[5] B.D. Jana and M.M. Stack, Wear 270 (2011) p.439. 

[6] Y. Petrov, V. Bratov, G. Volkov and E. Dolmatov, Advances in Rock Dynamics and Applications, 
pp. 163-183, Y. Zhou and J. Zhao Eds., CRC Press, Boca Raton 2011. 



11 



[7] J.R. Brockenbrough, S. Suresh and J. Duffy, Phil. Mag. A 58 (1988) p.619. 

[8] Y. Pctrov and N. Morozov, ASME J. Appl. Mech. 61 (1994) p.710. 

[9] N. Morozov and Y. Petrov, Dynamics of Fracture, Springer- Velrag, Berlin 2000. 



[10 

[11 
[12 

[13 

[14 
[15 
[16 

[17 

[is; 

[19 

[20 
[21 
[22 
[23 
[24 

[25 
[26 
[27 



Y.V. Petrov, Doklady Phys. 49 (2004) p.246. 

M.S. Binglcy and D.J. O'Flynn, Wear 258 (2005) p.511. 

A.P. Harsha and D.K. Bhaskar, Mater. & Design 29 (2008) p.1745. 

N.A. Kilchevsky, The theory of solid bodies collision, Naukova Dumka, Kiev 1969 (in Russian). 

I.Ya. Shtaerman, Doklady Akad Nauk SSSR 25 (1939) p.360 (in Russian). 

L.A. Galin, J. Appl. Math. Mech. (PMM) 10 (1946) p.425 (in Russian). 

F.M. Borodich and L.M. Keer, Int. J. Solids Struct. 41 (2004) p.2479. 

K.L. Johnson, Contact Mechanics, Cambridge Univ. Press, Cambridge, UK 1985. 

F.M. Borodich, Strength Mater. 21 (1989) p.1668. 

J. Woirgard, V. Auduricr and C. Tromas, Phil. Mag. 88 (2008) p.1511. 

Y.V. Petrov, N.F. Morozov and V.I. Smirnov, Fatigue Fract. Eng. Mater. Struct. 26 (2003) p.363. 

F.R. Tulcr and B.M. Butcher, Int. J. Fract. Mech. 4 (1968) p.431. 

Zh.-Ch. Ou, Zh.-P. Duan and F.-L. Huang, Int. J. Impact Eng. 37 (2010) p.942. 

R. Vignjevic, N. Djordjcvic, J. Campbell and V. Panov, Int. J. Impact Eng. 49 (2012) p. 61. 

G.A. Volkov, V.V. Silberschmidt, A.V. Mitrofanov, A.A. Gruzdkov, V.A. Bratov and Yu.V. Petrov, 
Doklady Phys. 54 (2009) p.322. 

LI. Argatov, N.N. Dmitriev, Yu.V. Pctrov and V.I. Smirnov, J. Friction and Wear 30 (2009) p.176. 

G.A. Volkov, N.A. Gorbushin and Yu.V. Petrov, Mechanics of Solids 47 (2012) p.493. 

V.I. Smirnov, Strength of Materials 39 (2007) p.46. 



12 



